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Abstract Using the hyperspherical adiabatic method 
with the realistic nuclear potentials Argonne V14, Ar- 
gonne V18, and Argonne V18 with the Urbana IX three- 
body potential, we calculate the adiabatic potentials 
and the triton bound state energies. We find that a 
discrete variable representation with the slow variable 
discretization method along the hyperradial degree of 
freedom results in energies consistent with the litera¬ 
ture. However, using a Laguerre basis results in missing 
energy, even when extrapolated to an infinite number 
of basis functions and channels. We do not include the 
isospin T = 3/2 contribution in our analysis. 
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1 Introduction ical orthogonal basis is flexible enough to handle the 

complicated nonadiabatic couplings. 

The hyperspherical adiabatic representation is well-established 


in atomic physics m ■ Few-nucleon problems have also 
taken advantage of this technique in solving simple model 
nuclear potentials [3]. Though convergence with respect 
to the number of included channels is typically favor¬ 
able in short-range low-energy atomic calculations [4], 
we set out to check the convergence in the three-body 
nuclear problem. In particular, we test whether a typ- 
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The solution of the adiabatic Hamiltonian, eigen¬ 
values (adiabatic potentials) and eigenvectors (channel 
functions) enters in the construction of a set of hyper- 
radial equations that can be solved for bound or scat¬ 
tering states. This gives to the hyperspherical adiabatic 
representation a great flexibility. Most of the applica¬ 
tion of this technique has been done to determine bound 
state solutions, however in recent years several appli¬ 
cations to determine scattering observables has been 
studied as well [o . 

In the present work, we use a hyperspherical har¬ 
monic basis to diagonalize the adiabatic Hamiltonian 
using the Argonne V14 (AV14) [6], Argonne V18 (AV18) [7], 
and AV18 with the Urbana IX (UIX) [5 three-body 
nuclear potentials. As a first step, to solve the adia¬ 
batic Hamiltonian we use a basis of symmetrized hy- 
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perspherical harmonics. In a second step, we solve the 
set of coupled liyperradial equations to determine the 
three-nucleon bound state using either a Laguerre poly¬ 
nomial basis or discrete variable representation (DVR). 
We compare the effectiveness of the two bases. This 
work can be seen as a first step in the application of 
the DVR technique, currently used in atomic physics, 
to the three-nucleon problem having in mind, as a final 
step, the treatment of the continuum spectrum. 

Section [2] describes the hyperspherical coordinate 
description we use to solve the three-nucleon problem. 
Section [3] shows the solutions to the adiabatic Hamilto¬ 
nian. Section [I] compares two methods in determining 
the three-body bound state and Sec. [5] concludes. 

2 Theoretical background 

We concern ourselves only with the relative Hamilto¬ 
nian H re \. We recast H le i in hyperspherical coordinates [TJ 
[2] in terms of five hyperangles denoted by Cl and a single 
length, the hyperradius R. The relative Hamiltonian is 
then a sum of the hyperradial kinetic energy, the hyper- 
angular kinetic energy, and the interaction potential, 

= (1 > 

Here, m is the average nucleon mass and ft is Planck’s 
constant. The exact form of the square of the grand 
angular momentum operator A 2 depends on the choices 


of the Jacobi vectors and of the hyperangles and is not 
needed here. 

The solution Fe(R, Cl) to Eq. © is expanded in 
terms of the radial functions R~ 5 ^ 2 Fe^(R) and the chan¬ 
nel functions $ V (R] Cl ), 

*e(R, n) = R~ 5/ 2 ^ f Ev {R)$ v (R\ «)■ (2) 

v 

The channel functions at a fixed hyperradius R form a 
complete orthonormal set over the hyperangles, 

J dCl$* v (R-,Cl)$ v ,{R-Cl)=5 vv ,, (3) 

and are the solutions to the adiabatic Hamiltonian S7), 

H ad (R, Cl)$ v {R- Cl) = U V (R)$ V {R-, Cl), (4) 

where 

H ^~ A2+ R l 5/ * + v U R ,n). ( 5 ) 

After applying Eq. © on the expansion Eq. © and 
projecting from the left onto the channel functions, the 
Schrodinger equation reads 

(-^J!e +u ’W- E ) F *'W (6) 

- E ( 2P ( R ) ^ ( R )) p EV (R) = 0. 

V 1 ' ' 

The hyperspherical Schrodinger equation Eq. © is solved 
in a two step procedure. First, H a ^{R, Cl) is solved para¬ 
metrically in R for the adiabatic potential curves U V (R). 

In a second step, the coupled set of one-dimensional 
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equations in R are solved. In Eq. ©> p„„' Q w' 
represent the coupling between channels, where 


p v AR) 



(7) 


and 


Qw’{R) 



d 2 ^ 
dR 2 


n 


( 8 ) 


the second Jacobi vector, and the total orbital angular 
momenta, respectively. Additionally, the spin S '12 and 
isospin T 12 along the first Jacobi vector, and the to¬ 
tal spin S and isospin T are fixed. Other restrictions 
include that L 12 + S '12 + T 12 must be odd for antisym- 
metrization and L and S couple to J = 1/2, the total 


The brackets indicate that the integrals are taken only 
over the hyperangle fi with the hyperradius R held 
fixed. 


angular momentum of the triton. In general, the sub¬ 
orbital angular momenta and spin are not good quan¬ 
tum numbers. Thus, the basis is not orthogonal and the 
adiabatic Hamiltonian, Eq. ©. is solved via a general¬ 
ized eigenvalue problem at each hyperradius. 


3 Adiabatic potential curves 

We expand the channel functions, <H> v (R-,rt), using a 
basis of hyperspherical harmonics 4>km{ J~2). Here, K 
labels the grand angular momentum quantum number 
and M labels the different degenerate states for a fixed 
K. The hyperspherical harmonics diagonalize the hy- 
perangular kinetic energy operator, where K 2 4>km(^) = 
I\ (K + A)(f>K m( fl)■ The main challenge comes from cal- 


at a fixed hyperradius, which has been worked out by 
the second author using the technique of Ref. ^S]. 

The hyperspherical harmonics are chosen to have 
certain symmetry properties. For example, subsets of 
basis functions are chosen to have orbital angular mo¬ 
menta L 12 , Ti 2 , 3 , and L. These correspond to the or- 


Table[l]shows an example of a set of basis functions. 
The different quantum numbers are listed in each col¬ 
umn with a given set of channels labeled 1 through 18. 
The number of basis functions comes from the range of 
grand angular quantum numbers K for each channel, 
where K varies from K m i n , to ensure linear indepen¬ 
dence, to ATmax in steps of two to preserve parity. The 
total number of basis functions represented here is 617, 
227 have L = 0, 185 have L = 1, and 205 have 
L = 2. This particular selection of channels and values 
of A' max describes the three-nucleon bound state with 
an accuracy of about 1 keV [5], 

The solid lines in Figure QJa) are the adiabatic po¬ 
tential curves generated from the basis in Table[l]for the 
AV14 nuclear potential. The full ground state potential 


dilating the potential matrix elements 


bital angular momentum along the first Jacobi vector, 


curve can be seen in panel (b). There are a number 
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Channel 

Z/12 

L 12,3 

L 

S 12 

T 12 

s 

T 

K 

1 *- min 

^max 

1 

0 

0 

0 

1 

0 

1/2 

1/2 

0 

120 

2 

0 

0 

0 

0 

1 

1/2 

1/2 

4 

120 

3 

2 

0 

2 

1 

0 

3/2 

1/2 

2 

120 

4 

0 

2 

2 

1 

0 

3/2 

1/2 

4 

90 

5 

2 

2 

0 

1 

0 

1/2 

1/2 

10 

90 

6 

2 

2 

2 

1 

0 

3/2 

1/2 

6 

90 

7 

2 

2 

1 

1 

0 

1/2 

1/2 

6 

90 

8 

2 

2 

1 

1 

0 

3/2 

1/2 

4 

90 

9 

1 

1 

0 

1 

1 

1/2 

1/2 

6 

60 

10 

1 

1 

1 

1 

1 

1/2 

1/2 

2 

60 

11 

1 

1 

1 

1 

1 

3/2 

1/2 

10 

60 

12 

1 

1 

2 

1 

1 

3/2 

1/2 

8 

60 

13 

1 

1 

0 

0 

0 

1/2 

1/2 

12 

60 

14 

1 

1 

1 

0 

0 

1/2 

1/2 

8 

60 

15 

2 

2 

0 

0 

1 

1/2 

1/2 

16 

40 

16 

2 

2 

1 

0 

1 

1/2 

1/2 

12 

40 

17 

1 

3 

2 

1 

1 

3/2 

1/2 

10 

40 

18 

3 

1 

2 

1 

1 

3/2 

1/2 

12 

40 


Table 1 Channel quantum numbers used in the hyperangular 
basis. Li 2 , ^ 12 , 3 , and L , indicate the orbital angular momen¬ 
tum of the (12) pair, that of the third particle relative to the 
(12) pair, and the total. S \2 and T 12 indicate the spin and 
isospin of the (12) pair, while S and T are the total spin and 
isospin. K m i n (A max ) is the minimum (maximum) grand hy¬ 
perangular quantum number. The given values are used to 
generate data in Figs. Ha), [2] and [3] 

of lines that asymptotically approach the zero energy 
threshold, indicating that these curves represent the 
fragmentation into three free particles. Two lines ap¬ 
proach a negative energy threshold of -2.22MeV, such 




R (fm) 

Fig. 1 (Color online) Adiabatic potential curves. Panel (a) 
shows the adiabatic potential curves for the AV14 potential 
using the basis set in Table H The dashed lines show the 
known asymptotic limits for the two lowest potential curves, 
approaching the deuteron binding energy of -2.22MeV. Panel 
(b) shows the ground state adiabatic potential curves for 
three different nuclear potentials. 

that these channels indicate the fragmentation into a 
deuteron and a free neutron. The two lowest poten¬ 
tial curves are not converged, but the known asymp¬ 
totic limits are indicated by dashed lines. In practice, 
the potentials are smoothly connected to their known 
asymptotic behavior. The lowest potential corresponds 
to an s-wave configuration between the deuteron and 
the neutron, where there is no angular momentum bar¬ 
rier. The second lowest potential corresponds to a d- 
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wave configuration, where there is a significant angular 
momentum barrier. 

It is known that the hyperspherical harmonic basis 
works well at small hyperradius, but requires a large 
number of basis functions to reach similar convergence 
at large hyperradius. We have seen that with the in¬ 
clusion of more basis functions, e.g. going to 1052 basis 
functions with 362 at L = 0, 335 at L = 1, and 355 at 
L = 2, that the lowest two potential curves in Fig. DJa) 
are indistinguishable from their asymptotic behavior on 
the scale shown. Specifically, A' max is taken to be 150 for 
channels 1 to 6, 140 for channels 7 and 8, 110 for chan¬ 
nels 9 to 14, and 90 for channels 15 to 18. Figure mb) 
uses this larger basis to show the lowest adiabatic po¬ 
tential curve for three different nuclear potentials. Solid, 
dashed, and dotted lines are for the AV14, AV18, and 
AV18 with UIX potentials, respectively. The difference 
between the different nuclear potentials is most clearly 
observed in the lowest adiabatic potential. The inclu¬ 
sion of the three-body potential, e.g., produces the low¬ 
est potential curve. On the scale shown in Fig. Oja), 
there is no visible difference between the channels of 
the excited states. The next section describes the cal¬ 
culation of the bound state supported by the lowest 
adiabatic potentials. 


4 Bound state energies 

With the adiabatic potentials U„(R), we solve the hy- 
perradial Schrodinger equation Eq. ([6]) to determine the 
bound state energy of the triton. As a first approach, we 
use a seven-point finite difference method to determine 
the P and Q matrix elements, together with a basis 
of Laguerre polynomials in the hyperradius R. There 
are two convergence criteria to consider. First, suffi¬ 
ciently many Laguerre polynomials must be used to 
describe the hyperradial wave functions. Second, suf¬ 
ficiently many channel functions must be included to 
reach convergence. 

Figure[2Da) shows a convergence study for the AV14 
potential as a function of the number of Laguerre basis 
functions Nl with the lowest 41 channel functions in¬ 
cluded. The hyperangular basis used is that described 
in Table[Q The solid (dashed) line is a strictly quadratic 
(linear) fit to the data from 1 /Nj J = 0.03 to 0.12 (0.03 
to 0.06). The two fits are used to estimate the energy as¬ 
suming an infinite set of Laguerre basis functions. These 
extrapolations are calculated at different numbers of in¬ 
cluded channels Whannei and shown in Fig. [U(b). The 
upper (lower) points are from the quadratic (linear) ex¬ 
trapolations. Solid and dashed lines are linear fits to the 
data points of Fig. mb). Extrapolating to an infinite 
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Fig. 2 (Color online) Extrapolations used to estimate the tri¬ 
ton energy for the AV14 potential using a Laguerre basis in 
the hyperradius. The upper panel shows the change in energy 
as a function of the number of Laguerre polynomials Nl for 
41 channel functions. The solid lines show strictly quadratic 
and linear fits to subsets of the data points. The lower panel 
shows the extrapolated values for different numbers of chan¬ 
nel functions JV c h a imei- The upper and lower points (dashed 
and solid lines) correspond to the quadratic and linear fits, re¬ 
spectively. The lines show linear extrapolations to the infinite 
channel limit. 

number of channels gives an estimate of —7.671 (2)MeV 
for the triton binding energy. 

This technique seems to miss about lOkeV of bind¬ 
ing energy when compared with other estimates 110] , or 
worse if estimating the energy by the last data point and 
not extrapolating to an infinite basis. We suggest that 



Fig. 3 (Color online) Extrapolations used to estimate the 
triton energy for the AV14 potential using a DVR basis in 
the hyperradius with the SVD method. Solid and dashed lines 
show strictly quadratic and linear fits to the data points from 
l/f^channel = [0,0.04]. 

the error comes from the finite-difference method to 
calculate the P and Q matrix elements. The couplings 
show complicated behavior at small R that may not be 
captured by this technique. We propose using a discrete 
variable representation mm in R together with the 
slow variable discretization method (SVD) |13j . Thus, 
the P and Q matrices are not needed since the hyper- 
radial solution is solved in sectors that are exactly di¬ 
agonalized. The only additional requirement is that the 
phase of each diagonalized sector must match between 
sector boundaries. 

Figure [3] shows the energy as a function of the num¬ 
ber of included channels Whannei- For this data, we 
use a total of 40 SVD sectors, with 20 ranging from 
R = 0.2fm to 3fm and 20 ranging from R = 3fm to 
50fm. Each SVD sector consists of a DVR basis of five 
points. The hyperangular basis used is that described 
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in Table |TJ We find that the data changes in the fifth 
significant digit when we increase the number of DVR 
points to 10. Further, we find that the results are ro¬ 
bust given changes in the number of sectors, so long as 
the potential minima is covered by enough points and 
the maximum range extends to the asymptotic limit 
at about 20fm. The solid (dashed) line of Fig. [3] is a 
strictly quadratic (linear) fit to the data including all 
points where 1/-/V c h a nnei < 0.04. 

Compared to the Laguerre basis, we achieve a no¬ 
ticeably lower energy for the same number of included 
channels. Said differently, if we do not extrapolate to 
an infinite number of channels, the DVR energies would 
produce a better estimate than the Laguerre basis. The 
DVR basis is not a variational method, however, thus 
we include more channels in our DVR analysis (the left¬ 
most data point includes the 70 lowest channels). Re¬ 
gardless, the method shows a clear convergence and we 
estimate the triton binding energy to be —7.685(l)MeV. 

The analysis is identical for the other nuclear po¬ 
tentials, with Fig.[3]being qualitatively similar for each. 
Using the larger hyperspherical basis described at the 
end of Sec. [U we estimate the triton binding energy to 
be —7.685(l)MeV, -7.619(2)MeV, and -8.475(2)MeV 
for the AV14, AV18, and AV18 with UIX potentials, 
respectively. These results are in complete agreement 
with those given in the benchmark of Ref. na showing 


the great flexibility of the SVD technique to treat the 
complexity of the adiabatic couplings. 

5 Conclusion and Outlook 

In this paper we calculated the adiabatic potentials 
curves for the three nucleon problem using realistic nu¬ 
clear potentials in order to calculate the three-nucleon 
bound state. In calculating the triton binding energy, 
we find that the SVD method with a DVR basis gives 
results consistent with the literature mi , while using 
a naive finite-difference method leads to inconsistent 
results. This suggests that the channel couplings are 
handled well by the DVR basis. 

Though the binding energy of three nucleons can 
be calculated by different techniques (see Ref. mi and 
references therein), this method may prove useful in cal¬ 
culating scattering data. As a preliminary study, we at¬ 
tempted to propagate the R matrix in the SVD method, 
but the hyperspherical harmonic basis used does not 
produce convergence at large values of the hyperradius. 
This led to difficulties as we could not accurately con¬ 
nect the potentials and coupling matrix elements to 
their asymptotic behavior at large hyperradius in or¬ 
der to switch to a traditional R matrix propagation 
method. To get around this limitation for scattering 
states below the three-nucleon breakup, we could use 
the method described in Ref. mi, which only requires 
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information from the bound state calculations in order 
to estimate the S matrix elements. For determining the 
breakup amplitude a larger basis should be necessary 
in conjunction with the method described in Ref. [5]. 
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